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Abstract: We determine a next-to- leading order result for the correlator of the shear 
stress operator in high-temperature Yang-Mills theory. The computation is performed via 
an ultraviolet expansion, valid in the limit of small distances or large momenta, and the 
result is used for writing operator product expansions for the Euclidean momentum and co- 
ordinate space correlators as well as for the Minkowskian spectral density. In addition, our 
results enable us to confirm and refine a shear sum rule originally derived by Romatschke, 
Son and Meyer. 
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1 Introduction 

The shear viscosity is one of the most important parameters characterizing the transport 
properties of the quark gluon plasma. Despite its demonstrated effect on the hydrodynamic 
expansion of the fireball produced in a heavy ion collision (see e.g. Refs. [1, 2] and references 
therein), its first principles determination has proved notoriously challenging. As lattice 
QCD is limited to the Euclidean formulation of the theory and weak coupling expansions of 
transport constants are tedious to perform and display slow convergence [3], much attention 
has lately shifted to gauge gravity calculations. This has in particular led to the conjecture 
of the quark gluon plasma being a nearly 'ideal' fluid, with a shear viscosity to entropy 
ratio close to the famous limit of 77/s > l/(47r) [4]. It would, nevertheless, clearly be of 
considerable value to have a quantitative estimate for the quantity starting from truly first 
principles, i.e. the QCD Lagrangian. 

At present, the arguably most promising approach to determine the shear and bulk 
viscosities in QCD utilizes lattice measurements of Euclidean energy momentum tensor 
correlators G n (x), which so far have been carried out only in pure Yang-Mills theory [5-7]. 
Specializing to the shear channel, n = rj, one uses this data in an attempt to invert the 
relation 

^r^> r^)H , 0<f<1> (L1) 

Jo 71 sinn 

to gather information on the spectral density, from the intercept of which the corresponding 
viscosity is obtained, rj = lim^^o Pv ■ To this end, Refs. [8, 9] derived a convenient sum 
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rule (cf. Eq. (1) in each reference and Eq. (4.13) below), which can be used to constrain the 
form of prj{w). This rule contains a contact term, which is obtainable from the UV-limit 
of a finite-temperature correlator that is most conveniently evaluated in the continuum, 
using perturbation theory and dimensional regularization. This calculation necessitates 
in particular the determination of the corresponding limit of the shear operator (T12) 
correlation function. 

In addition to sum rules, correlation functions of the energy momentum tensor are 
interesting in their own right. In Refs. [10-12], a systematic program was initiated for 
the next-to-leading order (NLO) perturbative evaluation of the Euclidean correlators and 
(Minkowskian) spectral densities corresponding to various components of the energy mo- 
mentum tensor in SU(iV) Yang-Mills theory. So far, the calculations have been limited to 
the bulk channel (the scalar and pseudoscalar operators), in which they have been applied 
first to the UV limit of the Euclidean correlator [11], and later to the full time-averaged 
spatial correlator [10] and the spectral density [12]. Direct comparisons with lattice data 
and gauge/gravity results (see e.g. Refs. [13, 14]) are possible in all of these cases, and are 
hoped to shed at least some light on the question, whether quark gluon plasma can be char- 
acterized as weakly or strongly coupled at temperatures slightly above the deconfinement 
one. 

In the present paper, our aim is to determine the shear correlator, defined in Eq. (2.10) 
below and referred to as the tensor channel correlator in some references, up to NLO 
in simultaneous weak coupling and UV expansions. The result of this calculation will 
enable us to verify the sum rule of Refs. [8, 9], in which we identify an additional order 
g 4 contribution that unfortunately is beyond the accuracy of our computation. We also 
evaluate the short distance limit of the equal time correlator in configuration space and the 
large frequency limit of the corresponding spectral function. The latter result will serve as 
a useful test for the prediction of Ref. [15], according to which the leading temperature- 
dependent corrections to the vacuum spectral density in the shear channel should vanish 
in the absence of fermions. Relevant work on the shear correlator has in addition been 
carried out in holographic models, cf. e.g. Refs. [16-18]. 

The outline of the paper is as follows. In Section 2, we explain our setup and define the 
quantities we set out to evaluate. Section 3 then contains a brief outline of the calculations, 
while in Section 4 we collect our results for the UV limit of the Euclidean correlator and 
apply them to the determination of the equal time spatial correlator as well as to the 
spectral density and the shear sum rule. In Section 5, we finally draw our conclusions. 
Some calculational details, in particular the detailed evaluation of one of the new sum- 
integrals encountered in the calculation, are left to the appendices. 

2 Setup and definitions 

Linear response theory describes the small departures from thermal equilibrium that arise 
when a physical system is subjected to an external perturbation. In particular, it relates 
changes in physical observables to the retarded correlators of the corresponding operators. 
When the energy and momenta are small, the behavior of the correlation functions can be 
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described using hydrodynamics, which can be viewed as a low energy effective theory for 
the underlying thermal field theory. The matching constants between these two theories 
are known as transport coefficients, and can be computed from the zero-momentum limits 
of retarded correlation functions using so-called Kubo formulae. 

The operator relevant for the computation of the bulk and shear viscosities is the 
energy-momentum tensor T^ v . The general tensorial structure of its correlators has been 
worked out in detail in Ref. [19], assuming symmetry in the indices, rotational invariance 
in D — 1 spatial dimensions, 1 and transversality with respect to the external momentum P 
(current conservation). However, depending on the definition, the correlator we compute 
may not be fully transversal. As pointed out by e.g. Refs. [8, 19], P^ (T^(P)T Q( g(— P)) 
may differ from zero by terms polynomial in P, which translate to contact terms in coor- 
dinate space. This is precisely what happens for the correlator conventionally computed 
in perturbation theory and on the lattice, and it is for this reason that we do not apply 
the general decomposition of Ref. [19] here. Instead, we project out the desired correlators 
without any reference to transversality. 

In this paper, we work with pure Yang-Mills theory at finite temperature. The dimen- 
sionally regularized Euclidean action is written as 

S E = J dr J d 3 ~ 2e x \F a ^F- u , (2.1) 

where we have defined F^ = (2i/g B )Tr {T a [D^, D„}} = 8^A a v - d„A* + g B f abc A^A c u and 
Dfj, = — ig B A^T a . The energy- momentum tensor of the theory reads correspondingly 

1 

and the correlators of this operator that we are interested in are defined through 



rp _ C rpa rpa rpa rpa try 



GWtfOc) = (T^(x) T a p(0)) c . (2.3) 

Here, the symbol (. . .) c stands for the connected part of the thermal correlation function. 
The corresponding momentum space expression reads 

G maP (P)= [ e- ip - x G^ aP {x) , (2.4) 

J X 

where the Fourier transform should be taken in D dimensions. 

To evaluate the shear channel correlator, Gi2,i2(^) 5 we introduce the projection oper- 
ator 

^HU,a/3 = 9 ^/MX^vfi + P^Pya) ' (2-5) 

where P^ u is a symmetric projector orthogonal to P as well as the four vector U = (1,0), 
P T = r P»P» + P 2 U^Uu ~ (UP) {UyP» + U»P,) (0 R , 



P P T — f) — IT P T P T P T — P T P T — D — 1 (2 7) 



x We work in D — 4 — 2e Euclidean space-time dimensions. 
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In components, Pq = = P^, = Sij — ^f. Applying X^ap to the correlator in 
Eq. (2.4) and choosing the spatial momentum p along the X£)-i-direction, we obtain 

X^af, G^ aP {P) = -D(D -2){D- 3) G 12>12 {P) , (2.8) 

where we have exploited rotation invariance in D — 2 dimensions. As suggested by this 
result, we define the Euclidean correlation function that we will evaluate below by 

G v (x) = 2c 2 v X^ a p(x) (T^x) T Q/3 (0)) C , (2.9) 

where X fJiV%a p{x) denotes the corresponding projector in coordinate space and we have 
(mainly for historical reasons; see the notation of Ref. [10]) introduced a coefficient c v , 
whose value can be chosen at will. When we set D = 4 and choose the spatial separation 
along the ^-direction, the correlator reduces to the simple expression 

G v (x) = -16cJ (T 12 (x) T 12 (0)) c . (2.10) 
3 Details of the calculation 

To evaluate the correlation function of Eq. (2.9) in momentum space, we use Eqs. (2.2)-(2.7) 
to obtain 

G V (P) = 2c% X ma0 [ e~ iP x (F- p (x)F- p (x) F^(0)F| (T (0)) c . (3.1) 

In expanding this expression to NLO in perturbation theory, we follow standard procedures: 
We start by generating the relevant one- and two- loop graphs, displayed in Fig. 1, with 
QGRAF [20], and next insert into them the Feynman rules and carry out the Lorentz and 
color algebra. Symmetrization and decoupling of momenta, i.e. removing as much of the 
numerator structure as possible, leaves now a considerably larger set of sum-integrals to 
be evaluated than in the bulk channel case. Keeping D unspecified for the moment, and 
using the gluon propagator in the covariant gauge with gauge parameter £, we obtain 

= D{D -f D ~ 3) (2 J? - J") -(D-2)(D- 3)JA + D(D - 3) (J* - J,\) 



, 2 N J D ( D 2 )( D 3 ) / 9J 000 , 4T 100 , 2T 010 , 4T 100 , i 27 h 

+ 9B iy c\ j I zx nooo + ^nooi + zx nooi + 4 - L i200i + iZX : 



010 
12001 



i 07-000 i 7-100 , /1 7-010 Q7-100 a 7-OOO a 7-IOO 1 7-OIO a 7-IOO 

T zx 11100 ' x 21100 "T ^21100 ~~ ^lOlOl ~~ ^lllOl ~~ ^11101 x 12101 ~~ ^111-11 

A 7-IOO 97-OIO 1 7-OOO 

4X 121-11 121-21 ' 11111 



D(D 2) / 200 o7 -020 , a 7-110 1 07-011 , 07-200 , 7-020 
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7-110 07-101 A (D — 2) 2 (£> - 3) / nr ooo n7 -ooo 07-020 

1X 12101 zx 12101 J 7 I lyx 12000 iyx 12001 ° x 12001 



+ D(D - 3) - 42^% + JOjffn + ^ I — D(D — 6) ( Zjtfn + 22 



-mil 



12 - WD + 3D 2 / 200 , ^ 002 \ Dp - 3)(3Z> - 10) 



— — 1 — I 07-200 , 7-002 \ , -v- ^, 100 ■ . . 
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(ii) (iii) (iv) 



O oo <X> 

(v) (vi) (vii) 



Figure 1. The LO and NLO Feynman graphs contributing to the correlators of the energy mo- 
mentum tensor. 

where the one- and two-loop sum-integrals J and X are given by 



[Pr(Q)] c [P 2 ] a+b -' 

q [Q 2 H(Q-Py 



Jab - 2P \riZ\a\trt D\21b ' \°-°) 



T f 9 h = j: [P T {Q)] S [Pt{R)] 9 [Pt{Q - R)^[P^+^+e-f- 9 -h-2 

abcde ~ Tq,r [Q 2 ] a [R 2 ] b [{Q - R) 2 ] C [(Q - P) 2 } d [(R - P) 2 Y ' { ' ' 

These integrals are defined using dimensional regularization in D = 4 — 2e dimensions, and 
we employ the shorthands P T (Q) = Q^Q V P%,{P) = q 2 - (q • p) 2 and f Q = T^ go J q with 
J q = f d jD ~ 1 q/(2-7r) D ~ 1 . Note also that the gauge parameter £ has cancelled in the sum of 
the diagrams, serving as a nice check of the computation. 

The evaluation of the specific sum-integrals present in Eq. (3.2) proceeds with methods 
closely analogous to those developed and explained in some length in Ref. [10] . One starts 
by performing the Matsubara sums, thus dividing the result into three parts proportional to 
zero, one and two Bose distribution functions. The next steps involve then an expansion of 
the result in positive powers of the momenta exponentially cut off by thermal distribution 
functions, and the subsequent evaluation of the remaining integrals. To highlight the 
differences with respect to the bulk channel calculation, we will in Appendix A go through 
the evaluation of one of the new sum- integrals, Zi2ioi> i R detail. For the others, we display 
only the result obtained after summing all of the terms in Eq. (3.2) together, cf. Eq. (4.1) 
below. This is in line with our general strategy of automatizing as much of the calculation 
as possible and all the time dealing with the full expression for the correlator rather than 
individual sum-integrals. 



4 Results 

4.1 The Euclidean momentum space correlator 

Following the steps outlined in the previous section, we arrive at a result for the sum of all 
integrals in Eq. (3.2). This leads to the Euclidean space expression for the shear correlator, 

G V (P) P 4 f _ 2 / A\ 2e (1 _ 13\ g 2 !^ (A\ 4e ( 4_ _ 206 24C(3) \ 1 
M A c 2 (4vr) 2 | 5 \P J \e 5 J (Air) 2 \P J \9e 135 5 J ] 
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9 (4tt) 2 



( 



22+ ^(t- p " 




2a7 [20 12 /p 2 2 



) 




qr 



™b(<?) n B (r) 



r 



) 



(4.1) 



in which p n is the Matsubara frequency and p = |p|. There are several interesting things 
to note in this expression. First, the T = part of the result agrees with the previously 
known expression obtained in Ref. [21]. Furthermore, as expected from the form of the 
shear operator, necessitating only coupling constant renormalization starting at NNLO, all 
1/e divergences have automatically canceled in the T-dependent part of the result. This 
in particular implies that there are no logs of P 2 in the thermal terms, which will have 
important implications on the spectral density, discussed in Section 4.3. 

We also note that the form of our result for the correlator allows an interpretation of 
its T-dependent part, denoted here by AG„(P), in terms of an operator product expansion 
(OPE). Following the notation of Ref. [10] (cf. this reference for additional discussion on 
OPEs), we obtain 



where bo = lliV c /(3(47r) 2 ) originates from the one-loop beta function of the theory, and 
we have used the thermodynamic identities 



There is an intriguing subtlety related to the renormalization scale running of the Wilson 
coefficients in the result. While the one multiplying the operator (e +p)(T) is expected to 
be scale invariant and thus contain logs of the renormalization scale A ~ oj at higher orders, 
the same does not apply to the trace anomaly contribution on the second row of Eq. (4.2). 
As discussed in Ref. [15], in order to obtain a transverse spectral function, this term should 
be interpreted as a genuine contact term, and the coupling in its Wilson coefficient allowed 
to run with a renormalization scale A, independent of u. This in particular implies that the 
NLO contribution to this function be a mere constant, which we have chosen to parametrize 
in terms of a coefficient £12 (cf. the constants and Q x defined in Ref. [10]). 

4.2 The equal time spatial correlator 

Next, we use the Euclidean momentum space result to study the corresponding equal 
time correlator in coordinate space, of which we are able to obtain the two leading terms 



AG„(P) 
4c 2 




(4.2) 




(4.3) 



(4.4) 
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Figure 2. The behavior of the Wilson coefficients 7 7); i(r') and 7r;;e+p(f) as functions of A^g-r, 
obtained using a two-loop expression for the gauge coupling g. The three curves correspond in each 
case to the variation of the renormalization scale A around the value 1/r by a factor of 2. The 
central value of the parameter is in accordance with the choices made in Ref. [10]; in our case, the 
method of fastest apparent convergence (FAC) is not applicable. 



in a short distance expansion. To this end, we take a four-dimensional (T = 0) Fourier 
transform of Eqs. (4.1)-(4.2), which using the standard formulae listed in Ref. [10] produces 

G v (r) _ 24d A . . , 2(e+p)(T) . . , (e-3p)(T) , , , / T 6 



4c^ 57rV 



with 



2(e+p)(T) (e-3p)(T) _ /T 6 \ 

+ 7,;e +P (r) + ^ lr,;e-s P (r) + O ^ J , 

(4.5) 



7,;i(0 = 1-?tS|+^ 4 )' ( 4 -6) 



9 (4vr) 

41 ff 2 jV, 
9 (4vr) 

7r?;e -3 P (r) = 0( 5 2 ) • (4 



7r);e+p (r) = 1 - - + , (4.7) 



The behavior of the Wilson coefficients 7n as functions of A^g-r is plotted in Fig. 2. 
Comparing the results to the corresponding bulk channel quantities from Ref. [10], we 
observe markedly improved convergence. It would be highly interesting to be able to 
compare the obtained behavior to lattice data, along the lines of Ref. [13] in the bulk 
case , as well as later for the full correlator, for which problems with lattice artifacts and 
discretization would be less pronounced. We are, however, unaware of the existence of such 
results in the shear channel. 
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4.3 The spectral density and the shear sum rule 

Another natural application of our results is the determination of the UV limit of the shear 
spectral density p n {oS) = ImG,(-i(w + ie),0), cf. the discussion in Ref. [10]. Analytically 
continuing Eq. (4.2) to Minkowski space and taking the imaginary part of the result at 
p = 0, we see that to the accuracy of our computation, no T-dependent terms survive, in 
accordance with the arguments of Ref. [15] . This leads to the compact result 

=-^^%;iH + (e+p)(T)% ;e+P H + (e-3p)(T)% ;e _3 P H + 0^) , 

(4.9) 



% ;1 H = l--^+0( ff 4 ), (4.10) 



where 

20 fN, 
9 (4vr 

7„ ;e+P M = 0( 9 4 ), (4.11) 
7„ ;e -3p(w) = 0(g 2 ) . (4.12) 

The computation of the full NLO spectral density is currently underway; to gain some 
insight into the features to be expected, we refer the reader to the bulk channel work of 
Ref. [12] as well as to the one- loop HTL result on Ref. [22]. 

The spectral density also appears inside a sum rule derived in Refs. [8, 9], which relates 
its integral to the energy density. In our notation, the rule reads 



1 f°° doj ( i 2 

-I^y_- W-)-^)l-o} + mn Q( U ,T) = - e(T) , (4.13) 

where the Minkowskian correlator in the contact term, Q(u),T), is defined through 2 

G(u,T) ee - AGv( ~^f = 0) -\J d 4 xe-o(r 00 (x)T o(0)) T , (4.14) 

and the subscript T denotes the finite-temperature part of the function. In Ref. [9], the 
entire contact term was argued to be finite and of order 0(1) x (e — 3p)(T). 

To inspect the sum rule, we begin by evaluating the contact term. To this end, we 
have independently checked the result quoted in Ref. [9], 

J d 4 xe^(r o(x)Too(0)) T = -jj{l + cV)}( e+ p)(T) (4.15) 

f l - g 2 b ln( Q0 )(e - 3p)(T) + o(g\ 1 



where Coo is yet another unknown three-loop parameter and where the renormalization 
scales are to be interpreted in the fashion discussed after Eq. (4.2). Adding this to the 



2 Note that when the corresponding Euclidean results are analytically continued to Minkowski space, we 
automatically take the real parts of the expressions. 
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UV limit of the shear correlator according to Eq. (4.14) and going to the oj — > oo limit, we 
obtain the simple (and finite) result 

hm Q(u;, r) = Iln^x(e- 3p)(T) = 0(g A ) , (4.16) 

which verifies the claim of Ref. [9] regarding the size of the term. The determination of 
the coefficients £12 and Coo is unfortunately beyond the accuracy of our calculation. 

Finally, we take a look at the integral containing the spectral function. For it, we may 
use the relation 

— \p v (u) - Pr ,{uj) | T =o ) = G„(0) - lim AG„(-iw, p = 0) , (4.17) 

in which we have used the fact that G(0) vanishes at T = 0. The first term here involves 
the zero momentum limit of the (Euclidean) shear correlator and thus reduces to two-loop 
vacuum type sum-integrals. This calculation is performed in Appendix B and results in 

which cancels with the leading order part of the contact term on the second row of Eq. (4.2). 
What remains from the sum of Eq. (4.17) is simply the (e + p)(T) part of the UV term, 
leaving us with 

= je{T) + 0(<, 4 ). (4.19) 

Here, the unknown 0(g 4 ) term contains contributions from both the In £12 term in Eq. (4.2) 
and the unknown three-loop contributions to the result of Eq. (4.18). 

Looking back to the left hand side of the sum rule of Eq. (4.13), we note that the sum 
of Eqs. (4.16) and (4.19) is in perfect agreement with the right hand side. As explained 
above, there however appear to be non-vanishing order <? 4 contributions to the former that 
our computation has unfortunately been unable to resolve. 

5 Discussion and conclusions 

In the paper at hand, we have perturbatively evaluated the UV limits of the thermal 
correlators of the shear operator T12 and the energy density Too i n hot Yang-Mills theory. 
The calculations were performed to next-to-leading order in both weak coupling and UV 
expansions, and were subsequently used to determine OPE coefficients for the equal time 
coordinate space correlator and the spectral density in the shear channel. In addition, 
we confirmed and refined a sum rule for the shear spectral density, proposed originally in 
Refs. [8, 9]. 

Let us briefly summarize our main results. In Section 4.2, we found that the NLO 
expressions we obtained for the Wilson coefficients of the coordinate space shear correlator, 
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Eqs. (4.6)-(4.8), display considerably improved convergence properties in comparison with 
those obtained in the bulk channel in Ref. [11]. For the spectral density, our Eqs. (4.11)— 
(4.12) on the other hand verified the claim of Ref. [15] that the leading order thermal 
corrections to the quantity should vanish for Nf = 0. And finally, in the context of the 
shear sum rule of Refs. [8, 9], 

~ 16^ /«, ^ " PlliU) |T= ° } = I ^ " <£%o SiU,T) ' (5 ' 1} 

we independently evaluated all of the terms, finding perfect agreement to the accuracy of 
our computation. In particular, we found a simple result for the contact term appearing 
here, cf. Eq. (4.16). The numerical evaluation of the constants £12 and £ooi defined in 
Eqs. (4.2) and (4.15), would unfortunately require a three-loop determination of the T12 
and Too correlators, which lies beyond the accuracy of our computation. 

There are clearly several directions, to which our present calculations can and should 
be continued. Most importantly, it would be highly interesting to perturbatively evaluate 
the full shear spectral density, using methods developed in Ref. [12]. This would provide 
important input for the lattice determination of the shear viscosity, as it would enable the 
analytic subtraction of the leading short distance terms from the lattice data and thus 
allow for a better modeling of the spectral shape. In addition to this, it would of course 
be important to consider the effects of fermions, both in lattice simulations and within the 
corresponding perturbative analysis. In the perturbative setting, the addition of dynamical 
fermions is in principle straightforward — though somewhat tedious — to implement; as 
a first step, it would be nice to be able to confirm the UV limit of the spectral density as 
predicted in Ref. [15]. And finally, one would naturally like to pursue the evaluation of 
the coefficients C12 and £00 , which are the only parts missing from the determination of the 
contact term in the shear sum rule. We plan to address these issues in future work. 
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A The two-loop sum-integral X^ 01 

In this Appendix, we present the detailed evaluation of the two-loop sum-integral 

which serves as an illustrative example of the subtleties related to having extra numerator 
structures in the integrand. Carrying out the Matsubara sums in this expression gives after 
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some straightforward algebra 



7-101 

x 12101 



P 1 



q,r Q 2 R\Q-R) 2 (R-P) 2 



Pt(Q)Pt(Q - R) 



+ 



n B (q) 



P 2 f q 



P 2 (l - r n /q n ) 



R \R 2 {Q - R) 2 (Q - P) 2 R 2 (Q - R)\Q - P) 2 
P 2 (l-p n /q n ) 



+ 



+ 



R 2 (Q - R) 2 (Q - Pf 
2P 2 

R 4 {Q-R) 2 (R-P) 2 

p2 



{Q-Py{Q-R) 2 (R-P) 2 



P T {R)P T {Q - R) 
Pt{Q)Pt{Q - R) 

P T (R)P T (Q - R) 



Q 



+ 



n B {q) n B (r) 



qr 



P 2 f q 



+ 



+ 



P 2 (l-r n /q n ) 
(Q - R) 2 {Q - P) 2 {Q-RY(Q-P) 2 

-Pn/qn) \ Pt{R)Pt{Q _ R) 

P t (R)Pt(Q + R) 
P t {Q)Pt{R) 



{Q-R)*{Q-P¥ 
2P 2 

{Q + R-P) 2 {Q-Py 

p2 



{R-Q-P) 2 {Q-RY 



(A.2) 



Q,R 



Here, the function /„ is associated with the appearance of squared propagators and reads 



fq ~ 2q* + 



1 + n B (q) 
2qT ' 



with 



J_d_ 

2q"dq 



n B {q) 



n B {q) 



(A.3) 



(A.4) 



\ E E {•••}• 

g„=±iq r n =±ir 

l)^ 1 and let Jq and J q stand for D and 



while the angular brackets stand for (cf. Ref. [10]), 

[■ ■ -]q = 2 E {•••}> [• • -\q, r 

q n =±iq 

As usual, we have here defined n B (q) = {e@ q 
D — 1 dimensional integrals, respectively. 

The parts of Eq. (A.2) proportional to 0, 1 and 2 Bose distribution functions n B (not 
counting those inside f q ) are henceforth referred to as the 0-, 1- and 2-cut contributions, 
and are dealt with using different methods. Of the three, the 0-cut piece corresponds to 
the vacuum (T = 0) correlator, and can be handled using standard integration-by-parts 
identities and integral tables, conveniently collected e.g. in the TARCER Mathematica 
package [23]. This calculation can thus be performed in a fully automated way. 

The 1-cut part, on the other hand, involves a three-dimensional thermal integral, 
in which the integrand is a one-loop vacuum amplitude. A crucial simplification in the 
evaluation of its UV limit, which we are presently interested in, is that we may perform an 
expansion in positive powers of the momentum Q, 

1 



(Q - R) 2 



J_ 2Q-R 4(Q • R) 2 
R? +_ W~ + ~R & + 



(A.5) 
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and similarly with R -H> P. This relies on the fact that Q is on shell inside the square 
brackets and the magnitude of the corresponding three-momentum is cut off by the Bose 
distribution function. The R integrals can then be performed using standard methods, 
taking advantage of identities such as 

P T (Q — R) = P T {Q) + P T (R) - 2Q„R v pT v , 

l(Q ■ P) 2 ] q = * 2 (^- 1 -pI), (q • p) 2 = ^ . (A.6) 

where we have made use of rotational symmetry. 

In the 2-cut part of the expression (A. 2), one may similarly expand the propagators 
V(Q ~~ P) 2 an d ~~ P) 2 i n inverse powers of P 2 . Defining z = q • r/(qr), q n = criq, 

r n = pir, with a, p = ±, we obtain 

1 1 



(Q - R) 2 2qr(pa - z) 
and subsequently end up dealing with the sums like 



(A.7) 



p=± a=± ' 

p=± <7=± ' 

1 y y - p = I y y ? — - = o , (a.io) 

4 ^-f ^-f (per - z) n 4 ^—f ^-f (pa — z) n v ; 

p=± a=± vl ' P=± <J=± w ' 

and their derivatives. The remaining task is then to perform the angular integrations. The 
terms odd in q and r vanish due to antisymmetry, and fixing the direction of r, we can 
write the remaining ones in terms of z- averages, a la 

,2 



q • r q 

1 - z 2 /- ~ 



(1) , (A.ll) 



where we have made use of rotational invariance and the dimensionally regularized angular 
integration measure. The other z-averages encountered in the calculation are 

1 \ D -3 , , , , , 

T^) z = D^ {lh (A ' 12) 

Xl-z 2 ) 2 / z - (D- 4)(D-6) (A ' 13) 
while the averages over the angle between r and p we need read 

<(* • p) 2 > ? = ^i(i>. <( f • p) 4 > f = a) • ( A - 14 ) 

Finally, for clarity of presentation, we wish to remove all n\(q) terms from the result, 
which is most conveniently done by utilizing the integration by parts identity 

nn{q)q n - l J q = D ~* + n [ n B (q)q n ~ 3 . (A.15) 
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Setting D = 4 — 2e, this leads us to the final result for the UV limit of the sum-integral 
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where we have in the last stage performed an expansion in e and switched to the MS 
scheme. Results for the other sum- integrals, as well as for the full correlator of Eq. (3.2), 
can be obtained in a fully analogous way. 

B The zero momentum limit of the correlator 



In this Appendix, we determine a perturbative result for the quantity G v (0), needed in the 
discussion of Sec. 4.3. In the zero momentum limit, the expression (3.2) consists of vacuum 
sum-integrals only (recall the definitions of Eqs. (3.3) and (3.4)) and collapses to 
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Using rotational invariance as well as symmetrizations, the integrals involving the trans- 
verse projection operator (i.e. those having non-zero upper indices) can furthermore be 
related to the more elementary quantities 
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21a ' abc 
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via the explicit relations 
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(B.4) 
(B.5) 
(B.6) 

(B.7) 

(B.8) 

(B.9) 

(B.10) 

(B.ll) 

(B.12) 



Using integration- by-parts (IBP) relations on the spatial momentum integrations in a sys- 
tematic way [24, 25], all two-loop vacuum sum-integrals 5~L reduce to products of one- loop 
ones I*, which by the same strategy (or by direct recursion = 2a ^~ Z) l h a following 



from the known analytic solution) can furthermore be related to a few one-loop master 
sum- integrals. Letting finally D — s> 4 and plugging the result into Eq. (B.l), we obtain 



G„(0) _ 8d A g 2 N c f n B (q)n B {r) 
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